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O ■ Abstract 

q-( In the two previous papers of this series we defined a new combinatorical approach to quantum 

gravity, Algebraic Quantum Gravity (AQG). We showed that AQG reproduces the correct infinitesimal 
dynamics in the semiclassical limit, provided one incorrectly substitutes the non - Abelean group 
SU{2) by the Abelean group U(l) 3 in the calculations. 

The mere reason why that substitution was performed at all is that in the non - Abelean case 
the volume operator, pivotal for the definition of the dynamics, is not diagonisable by analytical 
methods. This, in contrast to the Abelean case, so far prohibited semiclassical computations. 

In this paper we show why this unjustified substitution nevertheless reproduces the correct physical 
result: 

Namely, we introduce for the first time semiclassical perturbation theory within AQG (and LQG) 
which allows to compute expectation values of interesting operators such as the master constraint 
as a power series in h with error control. That is, in particular matrix elements of fractional powers 
of the volume operator can be computed with extremely high precision for sufficiently large power 
of H in the h expansion. 

With this new tool, the non - Abelean calculation, although technically more involved, is then 
exactly analogous to the Abelean calculation, thus justifying the Abelean analysis in retrospect. The 
results of this paper turn AQG into a calculational discipline. 
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1 Introduction 



In the two previous companion papers of this series PUI2] we introduced a new combinatorial approach 
to quantum gravity, called Algebraic Quantum Gravity (AQG) which uses ideas from Loop Quantum 
Gravity (LQG) [31 E]. One of the advantages of AQG over LQG is that semiclassical tools for 
background independent quantum field theories already available in the literature [I3 HE IE] can 
be applied also to operators encoding the quantum dynamics while in LQG this has so far been 
possible only for kinematical operators. The difficulty, as explained in detail in pQ has to do with the 
fact that in LQG the Hamiltonian or Master constraint operator [HI EH El Ej necessarily changes 
the number of degrees of freedom on which the semiclassical state, that it acts on, depends. The 
fluctuations of the degrees of freedom added by the operator are therefore not suppressed by the 
semiclassical state and their semiclassical behaviour is correspondingly bad. In AQG on the other 
hand the dynamics never changes the number of degrees of freedom and the just mentioned problem 
disappears. Furthermore, the annoying graph dependence of the normalisable coherent states of LQG 
disappears in AQG. 

In the companion paper [2] we have displayed a non trivial semiclassical calculation which shows 
that the Euclidean part of the Master Constraint for pure gravity in AQG reduces to its classical 
counter part in the H — > limit. The same calculation demonstrates that its fluctuations and quantum 
corrections are small. It is trivially extendable to arbitrary matter coupling and the Lorentzian 
constraint. However, the calculation was done using an, a priori physically unjustified, technical 
modification of the Master Constraint: namely we substituted the correct non - Abelean gauge 
group SU(2) of the canonical formulation of General Relativity by the incorrect Abelean gauge 
group U(l) 3 . The reason for why this was done is that the U(l) 3 analog of the volume operator 
|13j . without which none of the pieces of the Master Constraint (geometry and matter) can even be 
defined, can be diagonalised analytically. This is crucial in order that semiclassical calculations can 
be carried out. 

For the non - Abelean case on the other hand, except in special cases, the volume operator cannot 
be diagonalised analytically 1 , which prohibited so far any explicit calculations involving the quantum 
dynamics. This is one of the many criticisms that have been spelled out recently in |17j . 

In the present paper we show that this criticism, as many others, is unjustified, thereby intro- 
ducing semiclassical perturbation theory for AQG and LQG. Basically, what we are interested in are 
expectation values or matrix elements of (powers of) the Master Constraint Operator in semiclassical 
states. We will show how to compute those as a power series in h with error control! Notice that in 
usual perturbative QFT there is no error control. In fact, there the perturbation series is known to 
be only asymptotic for all realistic theories. 

As we will see, and as should be expected from [UI2], the non - trivial part of the corresponding 
h expansion consists in the evaluation of the matrix elements of fractional powers of the volume 
operator. The basic mathematical physics tools that we employ here are the spectral theorem for 
self - adjoint operators and the directed set structure of the cone of positive operators. More is 
not needed in order to define the series. The fact that makes the errors small on the other hand 
relies on the phantastically good semiclassical properties of the semiclassical states developed in 
For instance, if we are interested in cosmological questions, the power series is typically in terms of 
the classicality parameter t = AhG ~ 10~ 120 where A is the cosmological constant and G is New- 
ton's constant. This should be contrasted with the Feinstrukturkonstante a ~ 1/137 of QED. Here 
the extreme weakness of the gravitational interaction comes to help. Notice that our power expan- 

1 The matrix elements of its fourth power are known in closed form ^21 j however, the associated matrices on the 
invariant subspaces, while finite dimensional, are not diagonisable by quadratures. 
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sion, for cosmological settings, is as fast converging as for the new spin foam model introduced in [15J. 



This article is organised as follows: 



In section two we develop systematically semiclassical perturbation theory for AQG. For the 
benefit of the reader not interested in the detailed proofs, let us just state that if we are interested 
in the matrix elements < ip, V^ q ip' > of the volume operator V v associated with a given vertex v 
which is not calculable analytically, then we can replace V^ q within the matrix element, up to h k+1 
corrections, by the calculable quantity 



< Q v i> > 2q [i + E (-ir 1 9(1 q) f k q) ( - . ^— , - i)i (i.i) 

n=l 

where Q v is related to the volume by V v = \J~Qi- The operator Q v is a polynomial in flux operators 
and its matrix elements are known in closed form. 

In section three we show how the results of [H] can be used in order to perform concrete calculations 
(numbers!) using the tools of section two. This will show semiclassical perturbation theory at workl 

In section four we conclude and explain why the Abelean results of [2] carry over qualitatively 
and quantitatively to the non - Abelean case as far as the zeroth order in h are concerned while 
there are finite differences in the first and higher order corrections. This justifies the calculation of 
[2] and demonstrates that AQG is a calculational discipline. 



2 Semiclassical Perturbation Theory 
2.1 The Idea 

In semiclassical applications of AQG we are interested, in particular, in expectation values of (finite 
powers of) the Master Constraint with respect to coherent states ■0. As displayed explicitly in [T] the 
operators of which we have to take expectation values are then linear combinations of expressions of 
the form (symbolically) 



< ^, Pl {h)F 1 {V Vl )p 2 {h)F 2 {V V2 )..F N {V VN )p N+1 {h)tP > (2.1) 

where pj are certain polynomials in the holonomies along edges or loops adjacent to the vertices 
Vi, .., vn and Fj are certain functions of the volume operator V VI of the form Fj{V Vl ) = (Ql^ 1 ■ Here 
< qj = mj jni < 1/4 is a rational number and mj , rij are relative prime. The self - adjoint operator 
Q v is explicitly given by 2 

yj yj -yk yk yl yl 

jkl r> r> r> \ ' / 



where we have used the notation from PQE]: By ei a (v) we mean the edge of a cubic algebraic graph 
outgoing from vertex v into the positive (er = +) or negative (a = — ) direction respectively and Xj 
denotes the corresponding right invariant vector field of SU(2). Notice that ei + (v ) = ej(v), &i-{v) = 

2 For convenience we are working here with the dimensionless volume operator that results from the actual one by 
dividing it by a 3 , where a is the length scale parameter that enters the definition of the coherent states. 
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[ej(v — J)] x . The classicality parameter is t = l 2 P /a 2 where a is some length scale enters the definition 
of the coherent states. 



The naive idea to compute expectation values of (|2.2J) is as follows: Consider 



XI ■= 7 ,n VI i ^ ~ 1 (2 ' 3) 

The operator xj is bounded from below, xj > —1 and the quantity < ip,Q Vl ip > can be computed 
exactly by the methods of 5]. We have 

Fj(V Vl ) = | < V, Q„,V> > I 297 fiM, fj(xj) = (1 + (2.4) 

The idea is now to use the power expansion of the function t \— > f(t) = (1 + t) q , — 1 < t < oo given 
by 

q \ = +1 g(l-g)-(n-l + g) 

n / n! 



/W:=1 + E(0'" ( 

n=l v 7 v 



and to use the spectral theorem in order to get an expansion of Fj(V V[ ) in terms of xj of which 
coherent state matrix elements are computable by the methods of 5J, specifically 

/OO POO °° / \ 

mdEjit) = [i + J2 [I m ^/(t) 
1 n=l ^ ' 



n=l v 7 

where £7/ is the projection valued measure associated with xj. However, the second equality is wrong 
because the power expansion (|2.6J) does not converge outside the open interval t E (—1,1). Hence, 
this naive idea is false and must be substituted by a rigorous argument. 

Much of the effort that follows will be devoted to showing that the power expansion (|2.5|) can 
nevertheless be used in order to get reliable numerical results for the intended expectation value 
calculations including error estimates. More precisely, the power expansion up to order 2k + 1 gives 
reliable values for the matrix elements including h k corrections while the remainder can be estimated 
to be finite and of order h k+1 . The proof of this fact is somewhat involved and we therefore split it 
into several subsections. The reader not interested in the proof can jump directly to sectional 

2.2 Some Basic Tools 

In this subsection we will state and prove some basic lemmas which will be employed over and over 
again in the core part of the proof. 

Lemma 2.1. 

For each k > there exists < (3k < oo such that 

WW - ht 2k+2 < f(t) < Wt) (2.7) 
where fk{t) denotes the partial Taylor series of f(t) = (1 + t) q , < q < 1/4 up to order t k . 



4 



Proof. By Taylor's theorem we have 



/(f) = f k (t) + R h (t), f k (t) = 1 + ( I ) *"> Mt) = jf ds f{k+ ^ {s \ t - s) k (2.8) 
for any — 1 < f < oo. 

Since f^ 2k+2 \s) = -g(l - q)..(2k + 1 - g)(l + s ) q - 2k - 2 < for all s > -1 and (f - s) 2fc+1 is 
positive (negative) for t > s > (—1 < t < s < 0) it follows that R 2 k+i{t) < for all f > —1, hence 
f2k+i(t) > fit) for all f > — 1 and also for t = —1 by continuity. 

On the other hand, f( 2k+3 )( s ) = g(l - g)..(2fc + 2 - g)(l + s )9- 2fc - 3 > for all s > -1. Since 
(t — s) 2k+2 > for all s we certainly have f2k+2(t) < /(f) for f > and /2fc+2(f) > /(f) for 

-1 < f < 0. Now / 2fe+2 (f) = f k+1 (t) + ^ 2A , g + 2 ^ f fc+2 . Then for any (3 k > - ^ 2 ^ we still 

have / 2fc+ i(f) - P k t 2h+2 < /(f) for f > 0. Moreover, since < / 2 fc+i(f) — /(f) is bounded from above 
for -1 < f < there exists a finite (3 k such that also /2fe+i(f) - f3kt 2k+2 < /(f) for — 1 < f < 0. □ 

The optimal (lowest) value of is not easy to obtain by analytical methods but certainly by 
numerical methods. As an example, let us determine it analytically for the case k = 0. 

Lemma 2.2. 



l + gf-(l-g)f 2 < (l + f) 9 <l + gf (2.9) 
for all t > — 1 and all < q < 1. 

Proof. Consider /(f) = 1 + qt - (1 + f)". We have f(f) = g(l - (1 + f)^ 1 ) so that /' > for f > 
and /' < for —1 < f < 0. Thus, /(0) = is the absolute minimum of /. Let now g(t) = /(f) — j3t 2 . 
In order that g < we certainly need g(— 1) = 1— g — /3<Oso/5>l — g. We have 



</ = /'- 2/3* = ^t^p - 1) (2.10) 

where y — 1 + f > and p = 2/3/q > 2(1 — q)/q > 1 — g. Consider 

%) = ^ Zi T I =»/> 7 (y)= 2 _ q / lV1 fc(y), %) = (i-9)d/-l)-y 2 ^ + i/ (2.11) 
y - i y 2 q {y- i) 2 

It follows that k' = (2 — g)(l — y 1-9 ) hence j/ = 1 is the maximum of at which /c(l) = 0. Hence 
k < and thus /i' < is strictly monotonously decreasing. Since h(l) = 1—q and p > 1 — q it follows 
that g' is strictly monotonously decreasing for y > 1. The equation g' = or h = p has one more 
solution < < 1 because /i(0) = oo and h is decreasing. It follows that g' < for < y < y qtP 
and y > 1 and g' > for y q p < y < 1. Hence is a local minimum of g and y = 0, 1 are local 
maxima at which g(x = —1) = 1 — q — (5, g(x = 0) = 0. Thus, in order that g < it is necessary 
and sufficient that (3 > 1 — q where (3 = 1 — q is the sharpest bound. □ 

To see what lemma 12.11 is good for we state the following simple lemma. 

Lemma 2.3. 

Let B_ < B < E> + be self - adjoint operators and set B := [B + + B_]/2, AB := [B + - B_]/4. Then 
for any states ipi, i/j 2 in the common domain of all three operators we have 

|&(< ^i, [B - B]ip 2 >)|, |9f(< ij u [B - B\i) 2 >)| < < Vi, [A5]Vi > + < Vi, [AS]^i > (2.12) 
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Proof. Consider the the polarisation identity for any self adjoint operator B 

<^P 1 ,Btp 2 >=^J2 l< ^ B ^ > ( 2 - 13 ) 

(4=1 

where if) L = if)\ + «/V By assumption the operators B + — B, B — B_ are positive. Thus 

< g- > + l ~ < ^+ >-l ^ »/, , r? / . \ < g >1 ~ < g >-l . < ^+ >1 ~ < g- >-l 
£ < »(< ^1, >) = £ < 1 

< -B_ - < B + > +i ^. , „ , . < £? >_j - < B > +l < B + >_. ; - < B_ > +l 
I ± ±1 <r ox ^ „/, d„/,_ o _ i ^ + i +i 



| <3(<Vi,i^2>) = =^ ^< | 



;2.14) 



where < £> > t :=< if),,, Bip L > etc. We conclude 



|3£(<Vi,[a-% 2 >)| < < V'i, [Afi]^ > + < ^2, [A£# 2 > 
p(<^ 1 ,[ J B- J B]^ 2 >)| < <ip 1 ,[AB)ip 1 > + <if) 2 ,[AB)if) 2 > (2.15) 

as claimed. □ 

The idea behind this lemma is that, even if we cannot compute the matrix elements of B, provided 
we can compute the matrix elements of B and AB, formula (J2.15|) shows that the matrix element of 
B is a good approximation to that of B provided that the respective expectation values of AB are 
small. 

As an example, consider B = (Q 2 ) q with Q as in ()2.2|) polynomial in the elementary operators 
(flux operators for GR) so that matrix elements can be computed by the methods of [3] . We have for 
any A > that B = \ q [l + (Q 2 — A)/ A] 9 and obviously x = (X — A)/A > —1. Then we may apply the 
above lemma and set B + = A 9 [l+x], £>_ = \ q [l + qx — (1 — q)x 2 }. Given a coherent state if)\ we could 
set A =< ipi,Q 2 ip 1 > so that < ip h B+if)x >= X q and < ipi,B^ip! >= X q [l - (1 _ g )(<JkS£&> - 1)]. 
Hence, the expectation value of B = {Q 2 ) q is approximately given by < Q 2 > q up to a correction 
which is controlled by the fluctuation of Q 2 as one would expect. Notice that it has been established 
[S] that < Q 2 > agrees with the classical value to lowest order in h. 

2.3 The Expansion 

Turning to the general matrix element of interest ()2.1|) we expand 

< if), Pl {h)F 1 {v Vl ) P2 {h)F 2 {v V2 )..F N {v VN )p N+1 {h)if) > 

n 

= l[F°<if), Pl [l + h)p 2 {l + f 2 ]..[l + f N ]PN + Hp > 

n 

= HF^<iP,p lP2 ..p N+1 if)>+R} (2.16) 

where F® =< if), Q Vj ip > 2qj and the remainder R is a linear combination of terms of the form 

<if),p'J[..f l p' l+ ^> (2-17) 
with I = 1, ..,N and p'j are products of the P j while each /■ coincides with one of the fa. 
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Hence we are left with computing (J2.17j) . To do this, we use lemma 12.11 and find optimal /3j t k 
such that 

f~(t) := fok+iit) ~ P h kt 2{k+2) < W) < fi,2k + i(t) =: (t) (2.18) 
for all t > —1. We will now iterate lemma EPl frequently. In what follows we just consider the real 
part of ()2.17|) . the imaginary part works analogously as displayed in ()2.12|) . 

We start with 

< Tp,Plfl~fNPN+l1p > 
= < [P[(JV+l)/2]/[(iV+l)/2]-l- -flPl^}, /[(AT+l)/2][P[(AT+l)/2]-l/[(AT+l)/2]+l- JmPN+I^] > 
= ■ <lpl,f[(N+l)/2]lp2> (2.19) 

Here [t] denotes the Gauss bracket of the real number t. Application of ()2.15|) reveals that \R — Rx\ < 
R 2 + R3 where 

R = U(<1p,Pxf 1 ..f N p N+1 ll)>) 

R 1 = 3?(< Vij/kjv+ija]^ >) 
R 2 = < ipi, A/[(jv + i)/2]V>i > 

R 3 = <^ 2 ,Af [{N+1)/2] ^j 2 > (2.20) 

where fj = (fj' + fj~)/2 and Afj = (fj~ — /~)/4 are positive operators of which matrix elements can 
indeed be computed using the techniques of 0. 

Notice that in R\ we have achieved a reduction from N operator insertions /i,..,/jv of which 
matrix elements cannot be computed to N—l insertions /1, .., /[(at+i)/2]-1) f[(N+i)/2]+i, ■■, In- However, 
in Ri we still have iV — 2 (JV — 1) insertions if is even (odd) while in R 2 we still have iV (N — 1) 
insertions if iV is odd. It follows that by iterating ()2.20|) we can never achieve that all occurring 
expressions only contain operator insertions of which matrix elements are computable analytically. 
Notice that judiciously we started the iteration with the "middle" operator, otherwise we would 
create even more than iV operator insertions of which matrix elements cannot be computed. By 
choosing the middle, the number of those insertions at least does not increase. 

While the above procedure therefore never stops at a stage at which everything is exactly com- 
putable, we will show that after finitely many steps one reaches a stage at which the non computable 
terms can be estimated from above by computable terms and those estimates are of higher order in h 
than the order of the power expansion that we wanted to achieve. More precisely, we will show that 
after finitely many steps we obtain one expression which only involves the fj and thus the associated 
expectation value can be computed by the methods of |3] and contains h corrections up to order 
h k . In addition there are a finite number of terms each of which still contain at most iV operator 
insertions of which matrix elements cannot be computed plus at least 21 + 1 insertions of operators 
of the form A/j. These terms appear in the manifestly positive form, say for iV even 

<P{[AF[}P^..P([AF(]P 1 F 1 P 2 ..F N/2 P N/2+1 ^\ 

[af; +1 ] 1 p[[af[]p^.p;[af;]p 1 f 1 p 2 ..f n/2 p n/2+ ^ > (2.21) 

Here the F'-, F k are elements of the original set fi,..,f N . Likewise the Pj,Pk are elements of the 
original set pi, ..,pn+i- Instead of N/2 insertions of type F k we might also have some insertions of 
type Fk but since the latter have computable matrix elements, their estimate is of an even higher 
order in h than (j2.21|) as we will see so that ()2.21|) provides the type of term whose estimate gives 
the lowest power of h. Hence we will show that for I sufficiently large (J2.21|) can be estimated by a 
computable expression of order at least h k+1 . 
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2.4 Resolutions of Unity 

To estimate (|2.21j) we use the overcompleteness property of our coherent states [3] and insert resolu- 
tions of unity to cast (j2.21j) into the form 



du l}1 .. J dv x>N / 2 J dv 2 ,i.. J dv 2 ,N/2+i J dv' x x .. J dv' x l J dv' 21 .. J du' 2 l x 
x J du 3jl .. J d "3,N/2 J dv 4tl .. J du 4>N/2+1 J du' 31 .. J dv zl J dv' Al .. J dv' Al x 

JV/2 

X II < - 1 > < ^ 2 > fc ' -Pfc^l.fc > < lp2,N/2+l, Pn/2+\4> > X 

fe=l 

X <^2,l'[ AF /+lM,l > X 
Z-l 

x J] < Ik, [MjMj+i >< <j, Pjiij > < pWs,i >< tiu> i AF (}^ > x 

N/2 

X Ft < ^ 3 - fc ' F k^,k+y > < V>4,fc, -Pfc^3,fe > < 4>4,N/2+l, PN/2+l4> > (2.22) 
k=l 

Here the measures v^, are basically Liouville measures [S] over as many copies of the cotangent 
bundle T*(SU(2)) as there are edges involved in the expectation value under consideration and each 
of the integrals are over the entire corresponding phase space. The states ipi t k, ip'ij are labelled by 
the points in that phase space, see |H] for details. 

2.5 Sketch of the Estimate 

Before we give the rigorous argument, let us give its heuristic form in order to understand what we 
are driving at: In order to estimate (J2.22|) we make use of the following fact that has been proved in 

u 

< ij;i,Pi>2 >=< ^1,^2 > + hEkfaM] (2.23) 

where ipi,ip 2 are arbitrary coherent states whose overlap function < ipi,ip 2 > is sharply peaked 
(Gaussian) at the point in phase space by which they are labelled. The functions E , E\ are absolutely 
integrable against that Gaussian over both resolution measures u x , v 2 and are of zeroth order in h. 

Since (J2.23J) is sharply peaked at ipi = ip 2 we can basically set E (?pi,ip 2 ) = Sq^IjV'i) =< 
if>i, Pipi > in the essential support of the overlap function, the corrections being of higher order in h, 
just like the E% contribution. This is what results after integrating over the measure corresponding 
to ip 2 . Thus by integrating with respect to the measures u 2> k, v' 2 ji v ^k-, v' A j we can simplify ()2.22|) up 
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to h corrections to 



J du lt i.. J dv XtN / 2 J dv' l v . j dv' xl x 
x J dv^i.. J dv 3:N/2 J dv' 31 .. J dv' u x 



1-1 

xj] < i>[ d , [A^'W, J+ i >< % > < P'i^i >< lb, [AF/]^i,i > x 
i=i 



N/2-1 

x 1 [ < i>i,k, F k i/)i,k+i > < V'l.fc) -fV0i,fc > x 



fc=l 



X< 1pl,N/2,F N / 2 1p > < 1pl,N/2, PN/2^1,N/2 > < ^,P N /2+l^ > X 
X < [AF/ +1 ]^,! > X 

X J] < ^J, [A^^3J+l >< V>3,i, > < 1^1, P^l >< K^3,l > X 

J'=l 
iV/2-1 

X ] [ • lf)3,k,F k lf>3,k+l > < i>3,k,Pki>3,k > X 
k=l 

X < 1p3,N/2, FN/2i> > < 1p3,N/2, PN/2lp3,N/2 > < ^,Pn/2+1^ > 

In order to estimate (|2.24j) further we will prove that for arbitrary coherent states ipi,ip2 we have 
< Vi, [AF']^ >= h k+1 < ^,V 2 > [^(Vx,^) + fiG^Vi,^)] (2.24) 

where G' , G[ are absolutely integrable with respect to both ui, z/ 2 and of zeroth order in h. Similarly 
we will show that 

I < Vi,*V>2 > I < ^ 3 < >' [G (ip u ip 2 ) + hGxiMi)} (2.25) 

where Gq, G\ are both of zeroth order in h~ and both terms are separately absolutely integrable with 
respect to both V\, v 2 . Notice that while (|2.23|) and (|2.24j) are equalities, (|2.25j) is an inequality. It is 
due to this inequality, which provides the key part of the proof, that a negative power of fi appears in 
()2.25|) which is the price to pay for the estimate. The estimate is necessary to carry out because the 
left hand side is not calculable analytically. We conjecture that the matrix element (|2.25J) admits a 
sharper bound which does not involve the /i~ 3 power, however, the estimate presented below is not 
able to deliver such a result. 

Notice also the prime at the overlap function < ipi,ip 2 > in (|2.25|) . This prime is in order to 
indicate that < if>i,ip2 > is a Gaussian only with respect to the momentum variables of the phase 
space but not with respect to the position variables of six copies of SU (2) corresponding to the six 
edges adjacent to a vertex in the definition of the volume operator. This is unproblematic as far as 
integrability is concerned because the configuration space is a compact group and the Haar measure 
involved in v is normalised. However, the missing damping factor leads to an additional negative h 
power that comes from the measures v as we will see momentarily. 

Notice that the polynomials P are bounded operators so that | < if), Pip > | < ||P|| where ||P|| is 
the sup norm on the Abelean C*— algebra of functions of connections. Let us assume for simplicity 
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that ||P|| < 1 which is typically the case because in our application P is a matrix element of a 
holonomy in some representation of SU(2). Hence we can estimate ()2.24j) by 

dv^.. J dv 1)N / 2 J dv[ v . J dv'^x 

x J du 3tl .. J du 3tN/2 J dv' 3V . J dv' 3l x 
1-1 

x ft | < ^, [AF^' W > I I < ^, [AF/]^i,i > I x 

3=1 

N/2-1 

X I ) - 1pl,k, F k l(> 1>k+1 > | | < ^pl,N/2, FN/2lp > | X 
fc=l 

x| < [AF/ +1 ]^3,i > I x 

x J] I < Vsj, [A^14 +1 > I I < V4i, [Af^3,i > I x 

3=1 

AT/2-1 



X 



] [ < tp3,k, Fktp3,k+1 > | | < Tp3,N/2, F N/2 1p > 
k=l 

< ^(2i+l)(fe+l)-37V x 



x y dv 1A .. J du ltN/2 J dv' x v . J du[ t x 
x y dv 3>1 .. J du 3iN/2 J di/ 3 V . J di/ 3 l x 



l-l 

X 

3=1 

AT/2-1 

x 
fc=l 

X| < ^1,1^3,1 > ^ +1 (V> 3j l)| X 
l-l 

X II I < ^'^3,^+1 > ^(^+1)1 I < ^3,^3,1 > G' j(^3,l)| X 
3=1 
iV/2-1 



X 



II < ^3,jfc,^3,fc+l >' Gofc^fc+l)! I < 1p3,N/2, 4> >' G 0N / 2 (lp)\ (2.26) 
fc=l 

where Gofc(V') = Gok{ip,ip) etc. Basically, the remaining integral has as many Gaussians in the mo- 
mentum variables as there are measures. However, there are N factors of the form < ipi,ip2 >', hence 
there are 6N Gaussians in the configuration variables missing. Since the associated 6N resolution 
measures come with a factor of h~ 3 for each of which only fr^l 2 gets absorbed into the integration 
measure of the momentum Gaussian integral, the final integral in ()2.26|) is of order ft~ 3 / 2 ( 6Ar ) = h~ 9N . 
Hence the overall estimate is of order }i( 2l + 1 )( k + 1 )~ 12N a t most times a phase space point (defined by 
the state ip) depending factor whose value is of the same order of magnitude as the leading order of 
H coefficient of the calculable quantity 

< ^,F+..F+ 2 [AF{]..[AF; +1 ]^ > (2.27) 
10 



From this heuristic reasoning we conclude that in order that the correction terms be of order 
fjk+i we mus t have I > l]\f t k '■= QN/(k + 1). The actual value is slightly lower and is given by 
l Njk = (12N - 9)/(2(k + This holds for N > 2, for N = there is nothing to do and the case 
N — 1 can be reduced to the case N = by the expansion displayed in section El This means that 
the polarisation identity must be invoked at least l^,k times. The number of correction terms that 
are produced in each iteration triples the number of terms, hence the number of correction terms is of 
the order of 3 Zjv > fc and they are all of the same order of magnitude as the lowest order in h coefficient 
of (I2~2TI) . 

What all of this means is of course that up to a controllable error of order h k+l the value of (j2.19j) 
is given by 

< rl>,pifi.J NPN+1 ip >=< tl>,pitf..f+ PN+1 iP > +0(h k+1 ) (2.28) 

which is calculable and involves the fluctuations of Q up to order h k (we do not even have to calculate 
the constants In other words, we have shown that up to a controllable error of higher order in 
h it is allowed to use the power expansion of f(t) = (1 + t) q up to order 2k + 1 in order to get a 
reliable value of (j2.19|) . 



2.6 Rigorous Estimate 
2.6.1 Step I: Proof of (12331) 

We want to prove ()2.24j) and begin with the following observation 3 : 

Since (Q 2 - < Q > 2 ) 2 ( fc + 1 ) = (Q- < Q >) 2 ( fc + 1 )(Q+ < Q >) 2 ^ where < Q >=< ^j.Q^j > for our 
given coherent state ip we insert resolutions of unity for the matrix element of (Q— < Q >) 2 ( fc+1 ) 
between coherent states and can reduce our attention to the matrix element 

3 3 

< Vi, [Q~ < Q >]^2 >= ie jd2h []J < faitf'tii >~i[< ^ii,^2i >< >} (2-29) 

7=1 7=1 
= ie hh33 {[< ^11, y^21 > - < ^11, ^21 >< iflY^l >] < ^12, y 2 j2 ^22 >< ^13, *fV 2 3 > 

+ < 1pU,1p 21 >< ^y/Vl > [< ^12,^2^22 > - < ^12,^22 >< ^2^2 >] < ^13,^3^23 > 

+ < ^n,Vai >< ^iVi >< ^12,^22 >< faYi 2 ^ > x 

X [< ^13,^3^23 > - < ^13,^23 >< ^3^3 >] } 

Here ipi, ip 2 are any coherent states and the operators Yj for a given vertex v are defined by Yj(v) = 
t[Xj,(v) — Xj_(v)]/2, see (12. 2|) . where t = £p/a 2 is the classicality parameter for some length scale 
a. Here we have made use of the tensor product structure of the coherent states. It follows that 

< ip u , Y} ! i)' 2I > - < V17, ^27 >< faYpipi > 

= - { [< </>!,/+, X^;^ 2 ,7+ >~ <1pl+, Xj 1 + tfj I+ >< Vl,7+, ^2,7+ >] < 1pl,I-,1p2,I- > 

- < 1p 1>I+ ,1p 2 ,I+ > [< IplJ^Xfl^I- > ~ < •tyi- % X\ 1 jtyl- >< lk.,I-,lfa,I- >]} (2.30) 

We now have reduced the computation to matrix elements of right invariant vector fields on a single 
copy of SU{2) and can refer to [3] where we find 

it < ^ gi ,X^ g2 >= t < > & + — - coth(^ 12 )] Tr [f^ 2 ^ ) + 0(0 (2.31) 
£ Z12 sn^i 2 j 



3 The factor < , Qii>> is immaterial for the following power counting, because it can be factored out right fron the 
beginning. Therefore we will neglect it in the following discussion. 
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where z 12 is defined by ch(zi 2 ) := Ti{g^g')/2 for <7i,<?2 G SL(2, C). Thus 

lt[< 1P gi ,X 3 lPg 2 > - < 1P gi ,1Pg 2 X lj)g,X 3 ll)g > 



< Vfc, Vfc > [^) Tr ^^r,) - ^LjTttfgrj) + 0(t)} (2.32) 



where ch(z) := Tr(g j; g)/2. 

Thus, using resolutions of unity we have 



< V, (Q 2 - < Q > 2 ) 2 ( fc+1 V >= y ^i..^2( fc+ i) < A [Q- < Q >Hi > < fa, [Q- < Q >]^ 2 > •• x 
x < Vafc+i, [Q- < Q >]V»2(*+i) > < [Q+ < Q >] 2(fc+1 V > (2.33) 

Now consider for a = 1, .., 2/c + 1 with ip := ip 

\ dv a < V>a-1, [Q- < <2 >]^a > < Ipa, [Q~ < Q >}^a+l >= 77? TT / d ^Va,Ia I dfi H (U Ia ) X 
J t £ Jr 3 JSU(2) 

X eaq>(-[||p 0) j ff - Pa-l,I*\\ 2 + | \Pa,Ia ~ Pa+l,Ia\ ? + \\O a ,Ia - 9a-l,I*\\ 2 + \\0a,Ia ~ O a +l,Ia\ ?] A) X 
XPoU(z a<a - 1)Ia , Z a}a+ l tIa )Po\ 2 (z a ,a-l,Ia) ~ %Pla, Z a ,a+l,Ia ~ Ipia) (2.34) 

Here we have explicitly written out the resolution measures dv = cd?pd^n{U)\i + 0(h°°)]/t 3 where 
g = exp(pj(jj)U, U = exp(—i<Tj9 J ), c is a positive constant of order unity (essentially 7r~ 3 ) and have 
used the form of the overlap function < if) g ,i() g >= exp(— [| \p — p'\ | 2 + | \9 — 9'\ | 2 ]/i)[l + 0(ft°°)] close to 
g' = g where it is sharply peaked. Moreover, close to g = g' we have z = p+p' +i(9' — 9). The notation 
Pol n denotes a homogeneous polynomial of degree n. In order to evaluate (|2.34|) we use translation 
invariance of the Haar measure d^iH(U a ,ia) — dnH{Ua,,i<TU a -i,ia)~ l an d notice that close to 9 = the 
Haar measure d\in = sm 2 (9) sin(9)d9d(f)d(p with 9 2 = 9 2 , 9/9 = (sin(0) cos(^), sin(0) sm(ip), cos(0)) 
approaches the Lebesgue measure d 3 9. Let us introduce x 3 aIa := {p > aIa —p*)/s and y J aIa = (9 a j a — 
9 a -i,ia)/s where s = y/i. Notice that 9/s < ir/s but we can estimate the integral from above by 
extending the integral over 9 to infinity and in any case this just gives an error of higher order in h. 
Then we can perform the Gaussian integrals and find to leading order in t oc h 

/ du a < V> B -l, [Q~ < Q >]^a > < Ipa, [Q- < Q >}lpa+l >= t JJ X 
J la 

x exp(-3[||p a+ i )/CT -p a -i,ia\\ 2 + \\9 a +iM - 9 a ^ IcT \\ 2 } / {At)) X 

xPol^p, sx a _ hIa , sx a+1)Ia , sy a _ hIa , sy a+liIa ) Pol 2 (x a _ 1)Jcr , x a+1 j a , y a -i,ia, Va+ija) (2-35) 

where Pol^ is a not necessarily homogeneous polynomial of degree n. Hence after carrying out all 
2(k + 1) integrals we arrive to leading order in h at an expression of the form 

< ^i, (Q 2 - < Q > 2 ) 2{fc+1 V 2 >= t h+1 < Vi,V>2 > Poli 0(fc+ i)(p) Pol^ba,^,^,^) (2.36) 
as claimed. 

2.6.2 Step II: Proof of (|2T2"5|) 

In order to prove (|2.25|) we must estimate | < ipi, fip 2 > \- The method to estimate displayed below 
is not contained in j3] so that we must be more detailed here. In fact, the estimate (J2.25j) provides 
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the key part of the proof, all the manipulations with Gaussian integrals that follow later are fairly 
standard. 

Hence we need the explicit form of our coherent states |3] which for a single copy of 577(2) are 
given by 

= Ti^W & W = E (2-37) 

llwOII T 

where the sum is over all spin quantum numbers j = 0, 1/2, 1, .. while dj = 2 j + 1 and Xj is the 
character of the j— th irreducible representation ttj of 577(2). Relevant for a given function / of 
the volume operator located at a vertex u are only the states i/) gia with J = 1,2,3 and a = ± 
associated with the six outgoing edges adjacent to v. Here we use the notation of [TJ|2|, that is, 
the edges of the algebraic graphs are labelled as ej(v) =: ej + {y) denoting an edge outgoing from v 
into the I— th direction. The other three edges are ingoing but we can define the outgoing edges 
ei-(v) = [ej(t> — I)] -1 . This has the consequence that 

$ ei _ (v) (A) := 4> ei{v -i)(A) = ^dje-M+VfcMv - I))] mn [7r,-(A(e 7 (?; - J))" 1 )]^ 

j 

= E - /)) T ] m „T jmri (A( e/ _(i;))) 

$ ei+{v) (A) := &, W (A) = Ed i e- t ^' +1 )[7r j (^(i;))] mn [7r J -(A(e / (7;))- 1 )] nm 

j 

= Ev / ^e-^^ +1 )[ 7 r J (( 7/ (^))] m nT jm „(A(e /+ (t;))) (2.38) 



where ft, i— > Tj mn (h) = y/2j + l7Tj mn (h) denotes a spin network function and we have used unitarity 
in the second step. Let gj + := gi(v), := [giiy — I)] T then we obtain 



< 



^lj^2 > JJ||^l,/a||||^2,/a| 



l.a 



^2 n[ 2 ^ +i ] e "* i/£r(iiCT+i) e 

Oi<t} /cr |mj <r |,|mJ. .|,|Ti/ £r |,|n^|<j'i tr 
X Yl 7l jl,ra Ia n I A9l,Ia)^j ltT m' I<7 n' I<7 (g2,Ia) X 



x < ®/[Tj /+)mj+)n/+ <S>Tj I _ !mi _ tni _], f ®i \Tj I+ ,m' I+ ,n' I+ ® Tj I _ jm ' i _ tn ' i _] > 

^n^ +i ] e * fc+i) e x 



{jl<r} 1(7 \rnia\,\m' ItT \<jio 

X Yl 7r jl,m' Icr m I A92,Ia9ll a ) X 
la 



x < ®i[\ji+,mi + > <8> >], / ®j > ® >] > (2-39) 

where in the first step it was crucially exploited that the volume operator preserves the mutually 
orthogonal subspaces of the Hilbert space labelled by fixed spin quantum numbers. In the second 
step we exploited that the spin network matrix element in the third line is non - vanishing only 
when rij a = n' Ia and in that case is actually independent of nj a . The calculation then equals that in 
an abstract spin system so that we dropped the dependence on nj CT of the matrix element, thereby 
replacing Tj mn by \jm >. The latter fact follows for instance from the detailed analysis in 
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Denote the factor in the fifth line of (|2.39|) by ^4{ m)m '} and the factor in the sixth line by B{ m rn i} 
then ()2.39|) can be estimated by using the Schwarz inequality 



< E n^ +1 ] e " y " c, " +1) J E i^-'i 2 J E 

{jla} Ia V { m ' m '} V { m > m '} 



e n[ 2 ^ +i ] e ™^° /CT+i) jn^^w^^U) >< 

0/ CT } iff V /<T 

'E < > ® |ii-,mj_ >],/®j [b'i+,m J+ >® b'j_ s mj_] > (2.40) 

{m} 

where in the second step we used the completeness relation of the states \jm > on the fixed j sub- 
space of the abstract spin system. 

We will now estimate the second square root in (|2.4()jl . Since certainly f(t) = (1 + t) q — 1 < 
1 + (1 + 1) — 1 = 1 + 1 for alH > — 1 we have / < Q 2 / < Q > 2 and therefore have 

E < ®i[\ji+, m i+ >®\ji-,m I _ >},f<S>i[\ji+,m I+ >®|j / _,m/_] >< E < g ^2 m i° > II 2 

{m} { m } 

(2-41) 

where we exploited that in the angular momentum representation \j, m > = \j, —m > and that we 
sum over all m. Now we just have to evaluate the norm squared in (|2.41j) by using elementary angular 
momentum calculus, we do not need to work in the recoupling basis [13 . 
Since the Yj , Yj mutually commute for I ^ J we have 

Q 2 = —t 6 e IJK 6MNpYjYjY^Y^ [ Y^Yp 
= -SU'Y^Y^YJY^ 

= E { 

+[y/y/] [y]y}\ [y]y}\ \y*y*\ 

+[y/y/] [y]y]\ [Y*Yk] - [y}y?\ [y]y]\ [y*y£]} (2.42) 

When we take the expectation value of ()2.42|) with respect to ®i a \jia-, m ia >, om y the first and 
fourth term survive because all other terms contain factors of the form YjYf, YfYj, YfYf, YfYj 
of which Yj is diagonal but Yj, Yj is a linear combination of raising and lowering operators. 

Since %Yj = i(X J I+ — Xj_)/2 is represented on the abstract spin system as Jj + — Jj_ where J\ a 
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are the usual angular momentum operators we can now evaluate (|2.42|) further 

\\Q ®Ia \ jla,m Ia > || 2 
t 6 

= j ^2 < m i+i m i-\[ J i] 2 + { J iT + J t J i + J 1 r Jf\m I+ ,m I _ > x 

I^J,K; J^K 

x < m J+ ,mj_\ - [J|] 2 - [Jjf + Jjjj + Jjj}\m J+ ,mj^ > [m K+ - m K -f 
+ J < ™/+>™/-|[#F - [^7] 2 + [^7> J?]|m J+ ,mi_ > x 

x < m J+ , mj_ | [ J|] 2 - [Jjf + [Jj", JJ] |m J+ , mj_ > [m K+ - m K ^} 2 

I^J,K; J+K 

x{{ji+{h+ + 1) + Ji-tii- + 1) - ™/+ - ^/-)0'J+0*J+ + 1) + jj-0'j- + 1) - mj+ - m 2 j_) 
-{m I+ +m I _)(m J+ + m J _)} (2.43) 

where = Jj ±iJj = jf + - jf_ is a linear combination of the usual ladder operators and we have 
used the algebra [J 3 Ia , Jjf ,] = i 5 I j5 <7t7 >e jkl J l I<7 . 

The second term in ()2.43|) vanishes when summing over {m} while the first can be evaluated 
using the relation 

^m 2 = ^( J + l)(2 J + l) (2.44) 

\m\<j 

resulting in the compact formula 

E iiq ®*, i^,m /ff > ii 2 = {\m\\{2 3I<T + 1)] [jiE^^ + ^ ( 2 - 45 ) 

{m} la la 

As we wish to apply the Poisson resummation formula in order to evaluate the estimate ()2.4()j) . 
we should estimate ()2.45|) in terms of the integers ni a = 2ji a + 1. We have 

( 2 -mn(2j Ia + 1)] [ns>^ + 

' } la la 

J<r la 

1 t 6 

< (3) 5 Yll^+^- + n ?- n/ +] 

where we used y/mm' < (m + m')/2 and m 2 + (m') 2 < (m + m') 2 for non negative integers m, m'. 
Next we evaluate the characters 

I t \ sH^iaZia) , n A -. 

XjjA9l,Ia92,Ia92,Ia9l,Ia) = ( 2 - 47 ) 

where ^ 

ch(z Ia ) := -TT(g 1Ja gl I(7 g 2 j a g 1 j a ) (2.48) 
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Since the argument of the trace in (|2.48|) is a positive definite, Hermitean, unimodular matrix, it is 
clear that Zi a is positive. Since (|2.47|) appears under a square root which would make the application 
of the Poisson resummation formula cumbersome, we use the elementary estimate 



v / 2sh(^) < 2ch(z/2) (2.49) 
which holds for all z > 0. Putting everything together we find 

2t 3 e 3 */ 2 1 ^ 



x 

01,02,03 



X II{[E 4 ffj e^- /4 ch(n /iCT/ , /iCTI /2)] [ e-^^ /4 chK_ CT/ z 7 ,_ CTI /2)]} 

I n/ lCTj =l n/ i _ (T/ =l 

2 5 9V3VIl/ CT sh(^ 



,0-2,0-3 
00 



I ni >aj =— 00 n./ i _ (TJ -=-oo 

< t 3 e 3 */ 2 1 x ^ 

2 5 9v^vn/X^) Jb 3 x 

00 00 

X II^ S n^e-^VV'^W 2 ] [ ^ e -K-^/ 4 e «i.-x^,-^/2]} (2.50) 

I nj lCTj =— 00 nj i _ (T/ =-oo 

where we have used 2mn < m 2 + n 2 in the second step. 
Notice that 

II^H* . e-E«^^ - (2 ' 51) 

where ch(yjj a ) := Ti(gj Ia gjj a ) /2, j = 1,2. Using the Poisson resummation formula we can now 
evaluate 



n 

,2+ 



^e-V 2 = (2.52) 

n 

We-^V"/ 2 = ^e z2 /( 4i )(^ 4 + ^ 2 S 2 + ^ 4 /16) (2.53) 
^ — ' s y 2 4 
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where s = \/i and we have neglected 0(t°°) terms. We obtain the final formula 



| < rfr, frk > I < /TT sHyi ^ iy2 ^ TT e [^-^-M.xJ/W x 

x e nl a4+ ^ a+ ^/ 16 ] ( 2 - 54 ) 

We want to show that f)2.54|) is sharply peaked at Hij a = H 2 j a where gjj a = Ujj a Hjj a is a 
polar decomposition into unitary and positive definite Hermitean matrices of the SL(2, (D) matrices. 
Parametrising H = exp(pj<jj) = ch(p) + ajpj sh(p)/p where p 2 = pjPj and <jj are the Pauli matrices 
we easily find yjj a = 2pjj a and 

j i 

Pi IoV2 1(7 

<Mzu) = cHyi,io)ch(y 2 ,ia) + ci a sh(y 1Ia )sh(y 2) ia), c Ia = — — (2.55) 

Pl,IaP2,I<j 

Since ci„ G [—1, 1] we find from basic hyperbolic identities 

\Vl,Ia - y2,Ia\ < Zlo < + 2/2,/a) (2.56) 

Hence the argument of the exponent in ()2.54|) satisfies the following relation 

zj a - 2yl Ia - 2yl Ia < -(y 1Ja - y 2 ^f < (2.57) 

where equality is reached if and only if cj a = 1 and yij ff = y 2 j a , that is, Hi Ia = H 2 j a as claimed. 
This property of (J2.54)) will be crucial for what follows, because otherwise in the resolutions of unity 
that will follow below, if we would not have the Gaussian decay just established, the corresponding 
integrals would blow up. 

2.6.3 Step III: Proof of 

Next we prove ()2.23|) . From [5] we recall that the measure v involved in the resolution of unity is 
given, up to a 0(t°°) correction by 

du(g) = ^d 3 pdfi H {U), g = exp{p j a j ) U (2.58) 

where c is a numerical constant of order unity. Hence the resolution measure is up to the important 
factor 1/t 3 essentially given by the Liouville measure on T*(SU(2)). Also recall t 5j that the overlap 
function is essentially (i.e. up to 0(t°°) corrections) given by 

,9 |z|sh(p)sh(p') 4A 2 +28 2 +2<5 2 

A 2 = p2 + (p / )2 _~2 /2 



ch(2p) = (l + c)ch 2 (p + p') + (l-c)ch 2 (p-p')-l, c ' >,l[ 



ppi 

a. ft 

ch(z) = ch{p) cos(9) + zsh(p) sin(0) cos(a), cos(a 



PA 



p9 

5 2 = p 2 -~s 2 + 2 -9 2 , ~z = ~s + ij) (2.59) 
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where the notation is as follows: We set g = exp(pj(jj)U, g' = exp(p'j(jj)U' and exp(j)j(Tj) exp(p'-aj) = 
exp(pj<jj)U" and finally U = U^U'iU")' 1 = exp(—i9j(jj). Here 9 2 = 9j9j and all angles /3, a, <j> lie 
in [0,7r]. All quantities are uniquely determined by the above formulas. Moreover one can show: 
1. 

A 2 > where equality is reached if and only if pj = p 1 - . 
2. 

5 2 > where equality is reached if and only if either a) <fi = 9, \s\ = p, | cos(a)| = lorb)s = p = 
or c) (p = 9 = 0, Ti. 

Since 5 2 > 0, in order that the exponent in the first line of (|2.59jl vanishes we must have at least 
A 2 = 9 2 = 0. Since A 2 = is equivalent with pj = p'j this means that U" = 1 and therefore 6 = 
means U = U', that is, g = g' altogether. Clearly then z = p = 2p, hence = 9 = so we are in 
case c) and thus 5 2 = 0. 

What all of this means is that the overlap function essentially takes the form of a Gaussian peaked 
at g = g' . The more complicated structure is due to the fact that we are dealing with a non - Abelean 
group rather than a vector space. Since the overlap function vanishes rapidly when g ^ g', let us 
expand (|2.59jl in terms of Apj = p'j—pj and A9j = 9j — 9j where U = exp(—i6j(jj), U' = exp(— iO'jCTj). 
Then it is not difficult to see that 



< 1pg,1pg> > | 2 = C 

<Vi,M>| < ^r 3 e-^ [J ^[llpj + +pU (2.60) 



[Ap] 2 + [A9] 2 



9^ 



where Ap 3 Ia = - . 



A further result that we need from [Hj is the following: 

The polynomials P are typically just (products) holonomy operators in the spin 1/2 representation 
for which holds 

<VV,W 9 '> = < VV>VV > {[^ B ch(z/2) + (9'Tj)AB ^^ sh(z/2)] (2.61) 

+tl ^m + {ch{zW _ gt(gVw coth Ms h(;/2) - M + 
2z zsh(z) 

Here g^g' = ch(z) — irjZj/ zsh(z). Due to the Gaussian prefactor < ip g , ip g i > the first term of order 
O(t ) contributes only at g = g' and there is given by U ab where g = HU, g' = H'U' are the polar 
decompositions of g,g' with H = exp(ipjTj), U = exp(9jTj). At g ^ g' this term as well as the one 
of order t 1 and t°° are exponentially bounded in p,p' which is suppressed by the Gaussian. More 
precisely, in the vicinity of g' = g where the main contribution to the integral that we are interested 
in comes from, the term of zeroth order (in t) in (J2.23|) becomes approximately 

< V„ W 9 ' >= e-W S -vir*i*i-*if\ exp(z( P ; - Pj ) Tj /2) exp((^ + 9)^/2) (2.62) 

as claimed. 
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2.6.4 Step IV: Performance of the Resolution Integrals 



The way to extract the leading h power of the correction terms is now as follows: Recalling (|2.22|) . 
the final integrals to be computed are of the form 



dv\. 



dv 



N/2 



dv[.. / dv'j. 



du 00 .. / du 0N / 2 



'N/2 



dvL 



di>\.. \ dvi I du[ 



dv[ 



di>\ 



'01 



'ON/2 I dv m .. 

< 1p, P^QO > I I < ^00, Ml > I I < Vl> ^2^01 > 

< ^0N/2~lfN/2^N/2 > I I < 4 ) N/2, Pn/2+1^0N/2 > 

< lp, PiV>oo > I I < ^00, Mi > I I < 4>i, Pitioi > 

< '^N/2-lfN/2i ) N/2 > I I < 4 ) N/2,PN/2+li ) 0N/2 > 

< [Am > i x 

< P^' ol > | | < PW ol > | | < [A// +1 ]ft > 



du 01 . 



dv'oi x 



< lj>Q\,M2 > I •• X 



< ^01, M2 > | X 



<VoiJA/^> 
<^ 1) [A/^> 



< ^4-^02 > 



X 



< </4, P 2>02 > I - X 



(2.63) 



As we have seen, factors of the form | < if>i, fip2 > | are to leading order in h essentially estimated by 
products of Gaussians in p^ia ~P2,ia tl mes a homogeneous polynomial in p[ I(T , Vi,ia of order twelve. 
On the other hand, factors of the form < ipi, [AF]ip 2 > could be written to leading order in h as 
Gaussians in pP 2 Ia — p\ Ia and 6 2 Ia — 6\ Ia times a homogeneous polynomial of order 5(2fc + 2) in z\ a 
times a homogeneous polynomial of order 2k + 2 in z\ a — pP Ia where for gu a ~ Q2,ia we have 



'La 



[P2,Ia + Pl,Ia] + i[02,Ia ~ dl,Ia] 



(2.64) 



as follows from the explicit formulas proved in j^J. Here p 3 Ia = i < tp,X J Ia ifj > is the expectation 
value of the right invariant vector fields with respect to the state ip under consideration. Finally, 
factors of the form < ipi, Pip 2 > could be written, to leading order in h as Gaussians in p 2 Ia — p[ Ia 
and 6 2 Ia — 6\ Ia times sums of products of matrix elements of the SL{2, (D) elements exp(— i(j? 2 Ia — 
/'l./Jr- 2)., X p(^ /f , • 0\, ln )r j: 2). 

We now carry out the various integrals to leading order in h. We will drop the indices I,cr,j in the 
following power counting argument, hence for instance dv a means an integral over all the p J aIa , 3 aIa 

and by p a we mean the 18 variables p J a j a etc. and correspondingly ||p a || 2 := J2iaj(Pa,icr) 2 ■ Further 
dp a '■= Yliaj dp J a j a , dfin(U a ) := Ylia dfJ>H{Ua,ia)- Also, for simplicity we assume that, as it happens 
in our application, all volume operators are with respect to same vertex v, otherwise we just have to 
introduce more notation which however does not change the argument. Also for simplicity we assume 
that we just have to consider holonomies along the edges ej a (v), for loops we would otherwise just 
have to use more resolution integrals but they do not change the h power. 

The first type of integral we consider is (for a = 0, .., N/2 — 1 with ip := ip) 

dV0a\ < 4>a,Pa+l1p0a > \ \< 4>0a, F a +l4>a+l > \ 



< 



1 



dpo a djj, H {U { 



0« i 



exp(-[||p0o -Pa\\ 2 + \\P0a - Pa+l\\ 2 + \\0Qa ~ #a|| 2 ]A) X 



X 



t 3 J t 18 

|Poli 2 (Poa,Pa+i) Poli(exp([p 0a -p a ]/2),exp(i[6 l oa + a ]/2)| 



(2.65) 
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where Pol n denotes a homogeneous polynomial of degree n of the variables indicated and we have 
only kept the leading order in h. Here exp([p 0a — p a ]/2) stands for the collection of group elements 
exp(a J (y 0a /(T - etc. 

Using translation invariance of the Haar measure, we get d/j,H(Uo a ) = dfj,H(Uo a U~ l ) and since the 
Gaussian in (|2.64j) receives its essential support at Uo a = U a the Haar measure can be replaced by 
the Lebesgue measure, that is, for small Qi the Haar measure duuiU) = sin 2 (#) sm((p) d9d(j)d<p with 
9/9 = (sin(</>) cos(<f), sin(0) sin(yj), cos(0)) approaches d 3 9 with integration domain a ball of radius 
7i. The polynomial involving exp((# a + $a)/2) can be estimated from above by a polynomial in 
exp((p a — Pa)/2) alone because matrix elements of 577(2) group elements are bounded by 1. 

Let us introduce x 0a ■= {poa ~p)/s, x a := (p a -p)/s, yo a : = (9 0a - 9 a )/s where s = \/i then 
(J2.64j) can be further estimated by 



dv 0a \ < 1p a , P a +l4>0a > I I < ?Poa, fa+l^Pa+l > \ 

-72il dx ° a I dy 0a exp(-[\\xoa - x a \\ 2 + \\x 0a - x a+1 \\ 2 + ||yoa|| 2 ]) x 

1 7r 18 ^||j»a||<W 

x |Poli 2 (p + sx 0a ,P + sx a+1 ) Poli(exp(s[a;oa - x a )), exp(9 0a + 9 a ))\ 
< — / dx 0a / (i2/ aexp(-[||x a - a; a || 2 + ||x 0a - x a+ i|| 2 + ||?/oa|| 2 ]) x 

x|Pol 12 (p + sx 0a ,p + sx a+1 ) Poli(exp(s[xoa -x a )))\ 
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-Poli 2 (p)exp(-3(||x« - x a+1 \\ + s)74) (2.66) 

where we have used basic properties of Gaussian integrals, used f*' 8 dOI < J °° d9I for positive 
integrand J and dropped subleading orders of s. We will also drop the factor exp(— 3s||x a — x a +i||) 
in (J2.66|) in what follows as it just leads to higher order corrections. 

We can now perform the integrals corresponding to the measures v%, .., vn/2-i- Since the integrand 
no longer depends on Ui, .., t7v/2-i the Haar measure part of those measures just integrates to one 
and we are left with an integral of the form 

t 3N/2 g i8(A72-i) Pol6Af (P) J dxi '- J dx N/2~i exp(-3 1 \x - x x | | 2 /4) •• exp(-3 1 \x N/2 -i - x N/2 1 1 2 ) 
d 1 

= t 3JV/ 2 s i8(JV/2-i) Pq Mp) exp(-3| \x - x N/2 \\ 2 / (2N)) (2.67) 

where c, d are numerical constants and we have made use of the Markov property J dyK s (x, y)K t (y, z) ■ 
cK s+t (x,z) with Kt(x, y) = exp(— (x — y) 2 /t. The negative power of s = \/t comes from the fact that 
dv = dpdfi H (U)/t 18 but only 1/s 18 was absorbed by the Gaussian integral in p. 
Next we perform the integral corresponding to ^oat/2 which is of the form 

dpw/2d/iH(U N/2) 



dv N/2\ < ^N/2,PN/2+l%N/2 > I I < %N/2, [&F[Wl > I < C J £ig X 

x exp(-[||p 7v/2 -PN/2W 2 + ||poiv/2 -PiW 2 + [\\0on/2 - 9 N/2 \\ 2 + \ \9 0N / 2 - 9[\\ 2 ]/t) X 
x |Poli(exp((p iv/2 - Pjv/2)/2), exp(i(6» OA r/ 2 - 9 N/2 )/2) Poli ( fc +i)(Poiv/2 + p\ + i{Q[ - doN/2) x 
xPol 2(fe+ i)(p 7v/2 + Pi ~2p + i{9[ - 9 0N/2 )\ (2.68) 

We introduce new integration variables Xoat/2 = (PoN/2—p)/s, x[ = (p' 1 —p)/s, yoJV/2 = (6W/2 - #'i)/ s ! 
replace the Haar measure by the Lebesgue measure after an appropriate shift in SU(2) as above, 
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extend the integration domain from a ball to all of R 3 , estimate entries of SU(2) group elements 
from above by 1 and can estimate ()2.68|) further by 



dlS N/2\ < tpN/2, PN/2+lll>0N/2 > I I < *Pw/2, [^KWl > I < J dx QN / 2 dy 0N / 2 ) X 

x exp(-[||x 7v/ 2 - x N/2 \\ 2 + \ \x 0N /2 -x'i\\ 2 + [\\Von/2 + (0[ - N/2 )/s\\ 2 + 1 1 II 2 ]) x 
x |Poli(exp(s(x iv/2 - XAr /2 )/2)Polio(fc+i)(2p + s[x 0N/2 + x\ - iyoAr/2]) x 
xPol 2 (fc+i)(s[a;oiv/2 + x[ - iyw/2))\ 
= ct k+1 Po\ l0(k+1) (p) Po\' 2ik+1) (x[) exp([-3[||x' 1 - x N/2 \\ 2 + \\6[ - 9 N/2 \\ 2 /t]/A) (2.69) 

where we again used basic properties of Gaussian integrals and dropped subleading orders in h. We 
used the notation Pol^ for a generic inhomogeneous polynomial of degree n and performed estimates 
of the form |P 2 (fc+i)(z + w ) \ = P k+i{x 2 + y 2 ). 

Combining ()2.67|) and (J2.69J) we can perform the integral over v^/ 2 using the same manipulations 
and end up with an error estimate so far given by 

= ctfc+1 ~ 3AV2 g i8( N/2-1) PoW + io (fc+ i)(p)Pol 2(fc+1) K) exp(-3||x - x' 1 \\ 2 /(2N + 4)) (2.70) 
where x = p/s andp corresponds to the external state ip. Next consider the integral for (a = 1, .., I — 1) 

J dv'J < < , P^' 0a > I I < [A^a+i^+i > I (2-71) 
which can be estimated just like ()2.68|) with the obvious changes in the integration variables, resulting 
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f du'J < ^P'J' Qa > I I < [AF^' a+1 > I 

< ct fc+1 Pol 10 ( fc+ i)(p)Pol 2{fc+1) (x' a+1 ) exp([-3[||x' a+1 - x' a \\ 2 + \\9' a+1 - 6' a \\ 2 /t}/4) (2.72) 

We can now combine ()2.70|) and ()2.72|) and perform the integrals corresponding to the measures 
p'^ ..,u[_ v We use translation invariance of d[iH{U' a ) = ^/^([/^(L/^+i) -1 ) and introduce y' a = (9' a — 
Q' a +i)/ s an d arrive at the estimate so far 



^'^^^ ^^iv/^^ ^io'^M^tW)^) exp(-3||x - x[\\ 2 /{2N + 4/)) (2.73) 

where we absorbed the constants corresponding to Gaussian integrals over polynomials into c. 

All we did so far can also be done to the corresponding integrals involving the measures with a 
tilde resulting in the estimate 

^^^'^ ^s^^ PoW+^^ii^Pol,^,^^) exp(-3||a; - ~x' l \\ 2 /{2N + 4/)) (2.74) 

It remains to perform the integral 

1 



^(k+D-sN Pol 12N+wl(k+1) (p) J dv[J dv[J du' ol J du' ol x 

xPol 2(fc+ i)(^) exp(-3| \x - x' l \\ 2 /(2N + 4/)) Pol 2(fc+1) (59 exp(-3| \x - i' l \\ 2 /{2N + 4/)) x 

x I < P^' ol > I I < P^' ol > I I < f ol , [AF; +1 W ol > I (2.75) 
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We perform first the integral corresponding to ip' Ql using the same manipulations as in (|2.69|) leading 
to an estimate of ()2.75|) given by 

ct(2 m)( fc+ i)-3iv __L_ Po l 12jV+10(2m)(fc+1) (p) J du[j dv[J df>' ol x 

xPol 2(fc+1) (x;) exp(-3| \x - x' l \\ 2 /(2N + 41)) Pol' 2(fc+1) exp(-3| \x - x' l \\ 2 /(2N + 41)) x 

xPoi 2(fc+1) (x; +1 ) ex P ([-3[||x; - x ' 0l \\ 2 + - ^|| 2 A]/4) I < ^% > I 

< ^ (2 ^ +1)(fc+1) ~ 3JV ^8(4-1) P° W + 10(2Z + l)( fc+ l) (?) | d^/ X 

xPol' 2(fc+1) (x;) exp(-3| \x - x' l \\ 2 /(2N + 41)) Pol' 2(fc+1) (z|) exp(-3| \x - x' l \\ 2 /(2N + 41)) x 
xPol' 2{k+1) (x' l+1 ) exp([-3[||x;-^|| 2 + ||^-^|| 2 A]/4) exp([-[| \x\ - x' ol \ \ 2 + | \9[ - 9' 0l \ \ 2 /t]) x 
xPol 1 (exp((p; - p' 0l )/2), exp(i$ + 6' 0l }/2) (2.76) 

Using translation invariance dfiniU'i) = d^uiUliU^)^ 1 ) and c?/^h(C//) = ^^([//(f/^) -1 ) we can 
perform the Gaussians in 0| — 6^, fy' — 0q Z . After this, nothing depends on Q' ol any more and the 
corresponding Haar measure integrates to unity. Finally the remaining three momentum Gaussian 
integrals can also be performed and give a constant of order unity times s -18 from one missing 
Gaussian factor in the angle variables. Thus the final estimate is given by 

_ f (ai+l)(fe+l)-3iV-18(iV/2-l/2)p r . 1 / \ /q 77 \ 

The constant c is order unity because the resolution measures contan the right power of tt, see [3]. 
The bound depends on the point in phase space as expected. In order that the power of t equals at 
least k + 1 we should have I > (12N — 9)/ (2(k + 1)). Thus we see that we need, for given N the less 
terms in order to arrive at corrections of order h h+l the higher is the power of Q 2 — < Q > 2 in terms 
of which we expand the operator. 

Notice that the power counting argument outlined here can be supplemented by an analytical calcu- 
lation in order to obtain the actual value of the integral because after the estimate of < fip2 > 
everything is computable ^3]. However, the leading power in t of the integral will be at least the 
one we stated. Such an analytical calculation would also reveal all the finite, higher order (in h) 
corrections that we have dropped in the power counting argument. 

3 Explicit Example: The case N = 2 

Notice that for standard matter the terms appearing in the master constraint we have 2 < N < 12. 
We will exemplify the procedure for N = 2 and consider arbitrary k. This means that in order to 
obtain a computable error bound of order k + 1 we must have I > 15/(2(fc + 1)). 

We will work out only the real part, the imaginary part works completely analogously. We have 
with R =< iJ,pifip 2 f2P3^ > and ^i := Pi, 4>2 '•= ^2/2^3^ 

\®(R)-Ri\ < R2 + R3 (3.1) 

where 

R 1 = U{< ip,PifiP2f2P3ip >) 
R 2 = < V,Pi[A/]iPi^ > 

R 3 = < ^,P3/2P2[A/l]p 2 /2P3^ > (3.2) 
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The term R 2 is already computable by the methods of 5|. The first term involves the non - polyno- 
mial expression f 2 only linearly. Thus we can get rid of it by iterating the steps (j2.18|) and ([2. 19)1 . 
Hence we set ipi := p2fiPii ) , 1P2 '■= P3& an d find 

\R\ — -R4 1 < i?5 + Rq (3-3) 

where 

R A = U{< 1p,PlflP2f2P3lp >) 

R 5 = < i/j,pifip 2 [Af 2 ]p2fiPi'4> > 

R 6 := <^,p 3 [Af 2 ]p 3 ^> (3.4) 

All terms R4, R$, R§ are computable by the methods of [Hj. 

The term R 3 contains f 2 quadratically and iterating the steps (j2.18|) and (j2.19)l will not change 
that. However, we will iterate to generate terms of higher order. We set ip± = p 3 ip, ip 2 : = 
p 2 [Afi]p 2 f 2 p 3 %p and get 

\R 3 - R 7 \ < R 8 + R 9 (3.5) 

where 

R 7 = SR(< ll),p 3 f2p 2 [Af 1 ]p2f2P3ll) >) 

R 8 = <ij,p 3 [Af 2 }p 3 ^ > 

R 9 := <^,p 3 f 2 p 2 [Af 1 }p 2 [Af 2 }p 2 [Af l }p 2 f 2 p 3 ^> (3.6) 

The term R% is computable by the methods of [Sj while R7 involves f 2 only linearly and which thus can 
be gotten rid of by performing once more steps ()2.2()|) and (|2.21jl . Hence we set = p 2 [Afi]p 2 f 2 p 3 
and ip 2 = p 3 ip and find 

\Rr — Rio\ < R11 + R12 (3.7) 

where 

R 1Q = Sft(< 1p,p 3 f2p 2 [Af 1 ]p 2 f 2 p 3 1p >) 

R u = <4j,p 3 [Af 2 ]p 3 tp > 

R 12 := <^,p 3 f2P2[Af] 1 p 2 [Af 2 ]p 2 [Af 1 ]p 2 f 2 p 3 ij > (3.8) 

All terms R w , Ru, R\ 2 are computable by the methods of 0. 
Notice that we have shown so far that 

\R ~ -R4I < i?2 + -^5 + -^6 + -^8 + R\§ + R\\ + R\2 + R9 (3-9) 

where only the last term Rg is not computable. Since it contains three factors of the form Af 
meaning I — 1 in order that it already leads to order k + 1 corrections we must take k = 7. For lower 
values of k the procedure must be iterated again in order to generate more factors of the form Af 
in the non computable expressions. We will not do that as the general pattern should be clear by 
now. Notice that R 2 , R5, Re, Rs, -Rio, are of order h k+1 while Ru is of order ft 3 ( fc+1 ) and thus is 
subleading. 

We finish this section by noticing that in the special case N = 2 we actually do not need to 
use the iteration process in order to arrive at a calculable bound because we can use the Cauchy - 
Schwarz inequality to obtain 

I < *P,PlflP2f2P3^ > I < \\flPM\ \\P2f2P3ll>\\ 

< WhPnPW \\fmi>\\ 

< WffPMl II/2WII (3.10) 
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where in the second step we used that p 2 is a bounded operator with bound ||p 2 || < 1 an d in the 
third we used H/V'll 2 =< 4>, /V ><< i>, (f + ) 2 i> >< ll/^ll 2 - Both factors in the last line of ()3.10j) 
are computable by the methods of This will suffice to get a bound of order h but unfortunately 
not better than that. 

4 Conclusion 

In this paper we established that semiclassical calculations within AQG and LQG can be done 
analytically as a perturbation expansion in h. The errors can be controlled, i.e. they are finite and 
can be estimated from above and are of higher order in h than the order that one is interested in. 
This is despite the fact that the spectrum of the volume operator of AQG is not known analytically, 
hence disproving the negative claims made in "that nothing is calculable in LQG" . All one has 
to do in order compute expectation values, correct to the order h k , with respect to coherent states 
which involve the volume operator for a vertex v in the form {Ql) q , is to replace this operator within 
the matrix element by its power expansion up to order 2k + 1 that is 

< >* [i + 2 g (-D- rt'-y-«> (_!_, _ ir] (4 ,) 

The exact matrix elements of the operator Q v are known in closed form. The error term can be 
estimated from above and is of order h k+l . Hence we are in a situation similar to ordinary QFT, in 
fact the situation is even better because we have error control. 

Formula (j4.1j) is what one would naively do anyway in order to do semiclassical calculations, 
however, since the operator Q v is unbounded, one cannot use the spectral theorem to show that (j4.1|) 
really approximates the actual operator. In this paper we showed that the naive guess is nevertheless 
correct by using methods from functional analysis and properties of our coherent states. 

Notice that while we have done perturbation theory within AQG for cubic algebraic graphs, ev- 
erything goes through also in LQG for arbitrary graphs. Hence the present paper turns AQG and 
LQG into a calculational discipline. 
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